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Transport in photoactive graphene heterostructures, originating from the dynamics of photogen¬ 
erated hot carriers, is governed by the processes of thermionic emission, electron-lattice thermal 
imbalance and cooling. These processes give rise to interesting photoresponse effects, in particu¬ 
lar negative differential resistance (NDR) arising in the hot-carrier regime. The NDR effect stems 
from a strong dependence of electron-lattice cooling on the carrier density, which results in the 
carrier temperature dropping precipitously upon increasing bias. The ON-OFF switching between 
the NDR regime and the conventional cold emission regime, as well as the gate-controlled closed- 
circuit current that is present at zero bias voltage, can serve as signatures of hot-carrier dominated 
transport. 


Graphene, because of its unique characteristics, is of 
keen interest for optoelectronics research in areas such as 
photodetection, solar cells and light-emitting devices P- 
0 . Recently it was emphasized that graphene features an 
unusual kind of photoresponse mediated by photogener¬ 
ated hot carriers. The hot-carrier regime originates from 
the quenching of electron-lattice cooling when the system 
is close to charge-neutrality The hot carriers are 

exceptionally long-lived in graphene and can proliferate 
across the entire system 13©. This behavior, which sets 
graphene apart from other photoactive materials, leads 
to a dramatic enhancement in photo-response. 

The 2D character of electronic states, which are fully 
exposed in materials such as graphene, can enable new 
device architectures. One system of high current in¬ 
terest is stacked graphene-dielectric-graphene structures 
[see Fig§ a)]. Fabricated with atomic precision, such sys¬ 
tems can behave as field-effect transistors mm , reso¬ 
nant tunnel diodes (T2], [13], photodetectors [14j [15]; they 
also provide a platform to explore the Coulomb Drag ef¬ 
fect [16] and the metal-insulator transition D2I- In the 
hot-carrier regime, proliferation of photoexcited electron- 
hole pairs can result in an enhanced thermionic emission 
of hot carriers over the barrier. Vertical carrier extraction 
in such structures is facilitated by short interlayer trans¬ 
port lengths in the nanometer range, ultra-fast response 
times, and large active areas. Variability in properties of 
different 2D barrier materials (hBN, M 0 S 2 , WSe 2 , etc.) 
allows to tailor the photo-response characteristics to the 
different optoelectronic applications. 

Here we predict that interlayer transport in graphene 
heterostructures operating in the hot-carrier regime leads 
to an unusual type of photoresponse, namely, a neg¬ 
ative differential resistance (NDR). The mechanism of 
this NDR response relies on the interplay of two effects. 
First, the phase space available for phonon scattering 
rapidly increases with doping, enhancing the electron- 
lattice cooling and thereby altering the number of hot 
carriers in the system. Second, the large capacitance of 


the atomically thin device renders the carrier density in 
graphene layers sensitive to the interlayer potential dif¬ 
ference. These two effects combined together result in a 
reduction of the electronic temperature and a suppres¬ 
sion of thermionic current upon an increase of the bias 
potential V\>. The NDR effect arises when this suppres¬ 
sion overwhelms the increase in the field-effect transport 
under bias. Such an NDR mechanism manifests itself 
as an enhanced photo-current peaked at a bias potential 
well below the onset of the conventional field-emission 
regime, eVt> <C A, where A is the barrier height, as shown 
in FigfTJc). 

We note that the photoactive NDR architectures are a 
class of their own, and are well suited for optoelectronic 
applications. In particular, the fast response and in-situ 
tunability of graphene devices makes them ideal as photo¬ 
active switches or light-detectors with high gain. The 
NDR effect in photo-active devices, analyzed below, is 
distinct from the one in traditional NDR devices, such 
as Gunn diodes m or resonant tunneling diodes m 
[20] which rely on non-linearities under the application of 
large electric fields in the absence of light. 

Turning to the technical discussion, we note that, while 
in general both electrons and holes can contribute to 
thermionic transport, in practice transport is often dom¬ 
inated by a single carrier type. In the case of hBN, the 
barrier heights are A e i ^ 3.5 eV for electron transport 
and Ah ~ 1.3 eV for hole transport m ■ We can therefore 
treat the interlayer transport in an hBN-based system as 
dominated by a single carrier type (holes). 

Below we focus on the behavior in wide-barrier struc¬ 
tures, where thermionic emission of thermally activated 
carriers dominates over direct tunneling. This is the case 
for hBN thicknesses exceeding 4-5 monolayers (d ~ 1 nm) 
m at sufficiently high temperatures. Thermionic cur¬ 
rents are described by a particularly simple model when 
both graphene layers are at neutrality at = Q[22]: 

J(V b , T) = (<?o/e/?)e _/3A sinh(/?el/ b /2), (1) 
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FIG. 1. The Negative Differential Resistance (NDR) effect 
in a photoactive heterostructure operated in the thermionic 
emission regime. Shown here are (a) device schematics, (b) 
electronic band structure with the quantities discussed in the 
text marked, and (c) the I — V dependence under optical 
pumping obtained from Eqs. @- 0 . The bias voltage Vb 
controls the electron cooling through electrostatic doping of 
graphene layers. An enhancement in the cooling power upon 
increased Vb triggers carrier temperature dropping [marked 
T (l,2,3) 

in (c)]; see also simulation results in Figkla). The 
resulting suppression of thermionic current leads to negative 
dl/dV (in the grey region). For larger bias values, transport is 
dominated by field emission, yielding positive dl/dV outside 
the grey region. Shown here are results for both graphene 
layers undoped at Vb = 0. Results for nonzero doping are 
presented in Figs® 

where J is the current density per unit area, /? _1 = ApT 
and the go value is estimated in Eq. This expression 
follows from a general microscopic model at not too high 
bias Vb, such that the effect of anti-symmetric doping 
induced by Vb ^ 0 is stronger than the corresponding 
change in the barrier skewness [see derivation and dis¬ 
cussion in the paragraph before Eq. ©]• For larger bias 
values, field corrections to the barrier potential become 
important and must be accounted for; this is done in a 
microscopic model developed below. We also note that 
the electronic distribution is typically non-exponential 
when relaxation is slow. Slow relaxation would make the 
distribution tails more pronounced, ultimately enhancing 
the thermionic effects. 

The steep dependence of J on T and Vb in Eq. 0 
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FIG. 2. The closed-circuit current J cc induced by optical 
pumping (P = 10//W//im 2 ) in the absence of voltage bias, 
Vb = 0. In the transport regime dominated by hot carriers, 
the value and polarity of J cc are sensitive to the carrier den¬ 
sities rii, ri 2 in the graphene layers. Shown are experimental 
schematics (a), and the dependence J cc vs. n\ and U 2 (b). 
The four-fold pattern with multiple changes of the current 
polarity arises due to the strong cooling power dependence 
on carrier concentration (a non-linear color scale is used to 
amplify the features of interest). 

leads to NDR by the following mechanism. For electrons 
in thermal equilibrium with the lattice, Eq. 0 predicts 
a monotonic / — V dependence. A very different be¬ 
havior, which is key for NDR, arises under pumping. 
As pictured schematically in Figjljc) for three values 
T (b > T > T^ 3 ), in the hot-carrier regime the elec¬ 
tron temperature T becomes highly sensitive to Vb- The 
temperature-bias coupling arises because of the large in¬ 
terlayer capacitance producing bias-dependent doping in 
the graphene layers. An increase in carrier density leads 
to a faster electron-lattice cooling, which reduces thermal 
imbalance AT = T — To, with Tq the lattice temperature. 
The dependence in Eq. 0 then predicts suppression of 
thermionic emission. If strong enough, this suppression 
can lead to negative dl/dV. The NDR effect takes place 
in the grey region marked in Fig[l][c). 

The sensitivity of the electron-lattice cooling to car¬ 
rier concentration provides a smoking gun for the regime 
dominated by hot carriers, helping distinguish it from the 
conventional resonant tunneling NDR mechanisms, such 
as those discussed in Refs. m EG2- In that regard, we 
mention that for the more general case of unequal car¬ 
rier densities (n\ ^ ri 2 at 14 = 0)? our analysis predicts 
that the interlayer thermionic transport persists even at 
zero bias. As illustrated in Figj2j this produces a closed- 
circuit current at Vb = 0 with a characteristic density 
dependence: a four-fold pattern with multiple polarity 
changes. Such a pattern provides a characteristic signa¬ 
ture of hot-carrier dominated transport. 

The strength of the hot-carrier effects, reflected in the 
response of T to Vb [see FigQa)], can be characterized 
by the dimensionless quantity 

A dT 

a ~~Td(eV^y 


( 2 ) 
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The value a depends mainly on the power P pumped into 
the electronic system and on the thickness d of the barrier 
via the capacitance effect. As we will see, a governs the 
N-shaped I — V dependence: the condition for NDR can 
be stated as a > 1 / 2 . We will argue that values as large 
as a ~ 25 can be reached under realistic conditions. 

We now proceed to introduce our transport model. For 
thermionic transport over the barrier, as well as for tun¬ 
neling through it, we adopt a quasi-elastic but momen¬ 
tum non-conserving approximation. Indeed, a number 
of momentum scattering mechanisms at the interface are 
possible, such as scattering by defects, intrinsic phonons, 
substrate phonons, etc. Typical energy exchange in these 
processes is small on the barrier height scale A. With this 
in mind, we use the (quasi-elastic) WKB model for the 
interlayer transition matrix element: 

t(e) = r e -s /o'* P(*)y Ad = A - eE b x - e, (3) 

for e < A, and t(e) = T for e > A, where T is an energy- 
independent prefactor which depends on the barrier ma¬ 
terial properties. In Eq. ©, E\y is the electric field within 
the barrier due to interlayer bias, see Eq.([7|), m is the 
electron effective mass in the dielectric, and x* is the 
classical turning point for the skewed barrier potential, 
x* = min[d, (A — e)/eE^], see Fig. [ljb). Here, for the 
sake of simplicity, we ignore the effect of attraction to im¬ 
age charges, described by a —l/\x\ — l/\x — d\ potential. 
For a large barrier width, this gives rise to the barrier 
height Schottky dependence on the square root of E b. 
The effect is less dramatic for the not-so-large barrier 
widths analyzed below. 

Assuming an elastic but momentum non-conserving 
interlayer transport, the vertical current J{Vb) can 
be expressed through barrier transmission and carrier 
distribution [23]: 

J(V h ) = e^2 2 ^-\t ij (e)\ 2 D 1 (e)D 2 (e)[f 1 (e) - / 2 (?)], (4) 

e,ij 

where, due to the built-in field between layers, the en¬ 
ergy for the quantities in layer 2 is offset by e = e + eE^d 
[see Figjljb)]. Here J is the current per unit area, 
/i(e) = + 1 ] —1 , / 2 (e) = + l ] -1 are 

the Fermi distribution functions, the sum denotes inte¬ 
gration over e and summation over spins and valleys. 
We use tij(e ) = Sijt(e) defined in Eq.(|3|, and the den¬ 
sity of states per spin/valley Di^(e) = [e|/ 27 r (hv ) 2 with 
v ~ 10 6 m/s the carrier velocity. The temperatures es¬ 
tablished under pumping in unequally doped layers are 
generally distinct, (3i 7 ^ /? 2 , reflecting the cooling rates 
density dependence. Also, importantly, the electrostatic 
potential between layers, E\>d, is distinct from the bias 
voltage T4- This is so because the quantum capacitance 
effects, prominent at small carrier densities, [24] make the 
quantites /ii and /i 2 bias-dependent. These effects will be 
investigated below. 


A feature of our system which is key for NDR is the 
large mutual capacitance of graphene layers, which cou¬ 
ples Vb with the carrier density and makes the hot-carrier 
properties of each layer tunable. This coupling acts as a 
knob producing big changes in the hot-carrier photore¬ 
sponse through modest changes of carrier concentration 
on the order Sn ~ 10 12 cm -2 . Such bias-induced doping 
changes are routine in graphene/hBN systems [TO] . 

The bias-induced changes in carrier densities of 
graphene layers, as well as the electric field E b between 
the layers, can be described by a simple electrostatic 
model. We consider a dual-gated device with fixed charge 
densities in the top and bottom gates, respec¬ 

tively. The neutrality condition relates the charge densi¬ 
ties in the different regions of the device as 

n T + m + n 2 + n B = 0, (5) 

where ni, ri 2 are carrier densities on the graphene layers. 
The field E is related to these quantities through Gauss’ 
law: 

ftE h = 2ire (n T + n\ - n 2 — n B ), (6) 

where k is the dielectric constant of the barrier. A bias 
voltage kb applied between the graphene layers results in 

eV b = Mi - /J 2 + eE b d. (7) 

With doping-dependent fii,2> Eq.0 accounts for the 
quantum capacitance effects. Here we will use the T = 0 
expression ^ = sign(n*)/iu y^7r\ni\, which provides a good 
model over most of the relevant carrier density range. 

The three unknown variables ni, and E\> can now 
be found by solving the three equations , once the 

external variables tib and Vb are fixed. Through¬ 
out this work we focus on the symmetric case when 
tit = n b = —no with no interlayer bias applied (which 
corresponds to both graphene layers at neutrality when 
the gates are uncharged). In this case, Eqs. can be 

restated as n\ = no+£n, n 2 = no — Sn, and kE^ = AireSn. 
Then, plugging these values into Eq. ([ 7 ]), the density im¬ 
balance Sn can be obtained. The built-in field E\> matters 
in two different ways: the electrostatic potential value 
eE^d enters the WKB model, Eq. Q, as well as in the 
offset between D\{e) and ^(e) in Eq. ©>• 

The I — V dependence for the case no = 0 in shown 
in Fig[l] for finite values = n B = —no it is shown 

in Fi g|3} In our simulation, we first determine fi i ? 2 from 
Eqs.([5p([7| as a function of W and no- Using the (i\^ 
values, the electronic temperatures T\^ are determined 
from energy balance considerations, see Eq. (|9| below. Fi¬ 
nally, using the calculated values /iy 2 , Ti j2 , and E’b, the 
current is obtained from Eq. 0- We use hBN barrier pa¬ 
rameters for numerical estimates, with a (hole) barrier 
height A ~ 1.3 eV |2T, dielectric constant k ~ 5, thick¬ 
ness d = 6nm (~ 20 monolayers). Room temperature is 
assumed, Tq = 300 K, unless stated otherwise. 
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The prefactor value T in Eq. <© can be related to 
the measured conductance. Physically, T accounts for 
the processes in which tunneling couples to phonons or 
defects. The rates for these processes, which typically 
vary from interface to interface, can be estimated from 
transport measurements. Linearizing Eq. in eVb for 
ksT < A, and accounting for the thermionic contribu¬ 
tion due to e « A, the zero-bias conductance per unit 
area is 

G=^e~? A , g 0 = ^ND 1 (A)D 2 (A)\r\ 2 j, (8) 

where we take T\ — T^. Here N — 4 is the spin/valley 
degeneracy, and go is a prefactor in Eq. Q. The 
activation T dependence is consistent with that mea¬ 
sured for dark current |25]. Comparison with values 
G ~ 10 -7 Q~ 1 gm~ 2 measured at room temperature, and 
A ^ 0.4eV [25], yields values go ~ 1 Q~ 1 gm~ 2 and 
T~0.5eVA. 

Next we discuss how hot-carrier effects result in a cou¬ 
pling between the electronic temperature and the chem¬ 
ical potential for each graphene layer. In the contin¬ 
uous wave regime, the power P pumped into the elec¬ 
tronic system is distributed among the electron and lat¬ 
tice degrees of freedom. For simplicity, we will use a two- 
temperature model, describing electrons by a tempera¬ 
ture distinct from the lattice temperature, T > T 0 , valid 
when the carrier-carrier scattering rate is faster than the 
electron-lattice relaxation rate. Assuming spatially uni¬ 
form in-plane temperatures and chemical potentials, the 
total cooling power P obeys the energy balance condition 

P = P&c(m, Pi) + Popt(^ii Pi) + Pdiai^ii Pi): ( 9 ) 

written separately for each layer i = 1,2. Here we ignored 
effects such as direct interlayer energy transfer as well 
as the heat drained through the contacts. Equation © 
accounts for three cooling pathways intrinsic to graphene, 
mediated by acoustic and optical phonons (P ac , P 0 P t) [3 
6 , and the disorder-assisted acoustic phonon mechanism 
(Pdis) US , see Supplement. The cooling rates in Eq.([9| 
also depend on the lattice temperature To, however it 
suffices to treat To as a fixed parameter, since the heat 
capacity of the lattice greatly exceeds that of the electron 
system. 

The intralayer Joule heating is small and thus need 
not be included in Eq. ©• Indeed, typical vertical cur¬ 
rent values obtained in devices of active area ~ 1 gm 2 and 
under a bias Vo ~ 1V do not exceed a few nA mmm- 
This yields Joule heating sources which are at least three 
orders of magnitude smaller than the powers pumped op¬ 
tically. The interlayer energy transfer as well as the heat 
drained through contacts can be ignored in the energy 
balance in Eq. ([9| for similar reasons. 

The cooling power is a strong function of doping for 
the acoustic-phonon contributions P ac and Pdi S - This is 


J [nA/gm 2 ] 



FIG. 3. Map of the current J vs. eVb/A and zero-bias doping 
no (equal for both layers), for P = 10 gW / gm 2 . The current 
peak positions track the Vb values where the emitter is at 
neutrality (white dashed line). Local minima of J vs. Vo are 
marked with black dashed lines. The no = 0 slice, plotted 
in Figlllc), is marked with a dotted line. The peak-to-valley 
ratio (PVR) vs. no is shown in the inset. 

so because acoustic phonon scattering is dominated by 
quasi-elastic scattering processes at the Fermi surface, 
and also because of the strong dependence of the electron- 
phonon coupling on the phonon energy. The resulting 
dependence on the chemical potential takes the form[3 

nu 

P ac oc/z 4 (T-T 0 ), P dis (XAi 2 (T 3 -T 0 3 )/M (10) 

in the degenerate limit k^T g . The factor 1/kpi, 
where i is the disorder mean free path, describes the 
dependence of Pdi s on disorder strength. In contrast, the 
contribution P Qpt is essentially g- independent. Since the 
optical phonon energy in graphene is quite large, Tiuj o ^ 
0.2 eV, the value P op t is quite small, behaving as P op t oc 
exp(— huo/T) for &bT, g < Tuj o- The ratio P ac /P 0pt is 
therefore small near charge neutrality but can be order- 
one for strongly doped graphene with typical doping n ^ 
10 13 cm -2 [5]. The cooling power strong dependence on 
g can trigger the temperature dropping upon an increase 
in /i, i.e. dT/dg < 0, see Figga). 

The transition between the hot-carrier dominated 
regime, and the conventional field emission regime can 
be controlled by the power pumped into the electronic 
system. In our simulation we used the values for P typi¬ 
cal of laboratory lasers below saturation [27 (a few mW 
per a /im-wide spot). We take P to represent the power 
absorbed in each graphene layer (2.3% of incident power 
HD- The P values used are quoted in FigspFl panels 
and in Figs |2|3| captions. As shown in Fig§ c), a strictly 
monotonic I — V response obtained at P = 0, transforms 
into an N-shaped NDR dependence upon growing pump 
power. 
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The impact of doping no (taken to be equal for both 
layers at Vb = 0) on the I — V dependence is illustrated 
in Fig(3| The current peaks are shifted towards higher 
Vi > values upon no growing more negative [compare to 
Figg: c) which shows the no = 0 slice]. The peaks track 
the Vb values at which the emitter layer is bias-doped 
to charge neutrality. This is to be expected since the 
electron-lattice imbalance is maximal at neutrality. The 
peak-to-valley ratio (PVR) as high as ^ 5 can be ob¬ 
tained (see Fig. [3] inset). For no > 0, in contrast, the 
emitter layer is never at charge neutrality for any value 
of Vb, resulting in the N-shaped dependence fading out. 

The NDR effect is suppressed under a high bias poten¬ 
tial when field emission of carriers with energies below 
the barrier height overwhelms thermionic emission. As 
shown in Fig jTJ the high bias region eVb/A > 1 is charac¬ 
terized by J monotonically growing with increasing Vb. 
This behavior arises because lowering the barrier facili¬ 
tates tunneling and also because growing carrier density 
results in a faster cooling, thereby reducing the electron- 
lattice thermal imbalance [see Fig§ c) inset]. 

Next, we proceed to derive a simple criterion for NDR. 
We will focus on the fully-neutral case no = 0 (both 
graphene layers undoped at Vb = 0) pictured in Figjl] 
In this case, due to symmetry, we have T\ = T<i = T 
and ui = — ri 2 for any eVb- Further, assuming a small 
bias and/or a not-too-wide barrier, we can approximate 
the bias-induced chemical potentials as fi\^ ~ d=eVb/2. 
This simple relation is valid for eVb k(Tiv) 2 /e 2 d, cor¬ 
responding to the last term in Eq. § much smaller than 
g i — g 2- Lastly, accounting for the dominant role of 
thermionic emission, we model transmission as a step 
function, |£(e)| 2 « T 2 6(e — A). We integrate in Eq. Q 
over energies e > A max^T, Vb], approximating 
D i,2(e) as a constant and the Fermi distribution tail as 
e -(3(e-n)' This yields Eq. (|T|) to leading order in Vb/A 
and k^T/A with the prefactor go given in Eq.([8|. While 
the validity of Eq. 0 is limited to d which are not too 
small and also not too large, we find that it predicts NDR 
in the parameter range close to that found from the full 
microscopic model used to produce Figs§§ 

The criterion for NDR can be derived by taking the 
derivative dJ/dVb in Eq. 0 and setting it to zero, giving 

( 1 + ^) tanha; “/k = i’ X = W2 ’ (11) 

where a is the quantity — (A/eT)dT/dVb introduced 
above, describing the carrier temperature dependence vs. 
Vb- The a value controls the NDR effect. Maximizing 
the left-hand side in x we find the value /(A) = A 1 / 2 — 
(A —1) tanh -1 A -1 / 2 parameterized with A = l + (/3A) -1 , 
which is attained at x* = sinh -1 y/f3 A. It is straightfor¬ 
ward to check that /(A) < 1 for all A > 1. This gives the 
NDR condition a > 1/2. Below we use this condition, 
derived for no = 0, as an approximation for the more 
general case of no ^ 0. 


To estimate ccasa function of the model parameters, 
it is convenient to factorize a by applying the chain rule 
as a = qlt - ctn, giving 


olt{P) — — 


g dT 
T~dg’ 




A dg 
eg dV ]b 


( 12 ) 


Here olt depends only on the cooling pathways through 
Eq. |9|, while depends only on the barrier properties 
through the quantum capacitance effect of Eq. Q. Below 
we use Eq. (12) to estimate a and show that the NDR 


condition a > 1/2 can be readily met. 

To estimate olt we analyze the degenerate regime g 
k^T, where the doping-dependent contributions P ac and 
Pdis dominate over the roughly doping-independent P op t 
[see FigQa) inset]. In this regime, the cooling power 
behaves as Pi = 7 ig a (T b — Xq), with a = 4, b = 1 for 
acoustic phonon cooling, and a = 2, b = 3 for disorder- 
assisted cooling (here 7 i are constants that depend on 
the cooling pathways, see Supplement). We assume, for 
simplicity, that a single cooling pathway dominates over 
other pathways. Then Eq. (|9| yields 


a T = (0/6) [1 - (T 0 /T) b 


(13) 


This gives 0 < olt < a/6 with the low and high val¬ 
ues corresponding to T « To and T To, respectively. 
The crossover between these values occurs at a threshold 
pump power P* ^ 0.5 gW / gm 2 that marks the onset of 
the hot-carrier regime under typical experimental condi¬ 
tions. We define P* as the value for P at typical carrier 
densities n ~ 10 12 cm -2 , To = 300 K, and kp£ = 100 such 
that (T — To)/To = 0.1. This yields the above P* value. 
Maximum a T values found from Eq. (13) are 0.6 and 4 
for the Pdis and P ac pathways, respectively. 

Next, we estimate as a function of the barrier width. 
From Eq. 0, specializing to the case no = 0, we find 


~ 1. 


A/2 g 


6.5 


(4e 2 / nhv)kFd 1 + 0.3 d[ r 


(14) 


where fcp = g/Tiv. Here we have used the hBN barrier 
value A ~ 1.3 eV and g ~ 0.1 eV for typical bias-induced 
doping. This gives limiting values a /1 (d d*) ^6.5 and 

d*) = 0, with the crossover value d* ^ 20nm. 
From the above we see that the quantity a = ar • can 
reach values as high as a ~ 25. 

While the NDR criterion a > 1/2 is insensitive to 
the cooling mechanism (so long as it is a strong func¬ 
tion of carrier density, as discussed above), the form of 
the I — V dependence may reflect the cooling mecha¬ 
nism specifics. This is illustrated in FigQb) for the 
Pdis mechanism. In particular, we consider the bias E p k 
where the current peaks. From Eq. ( pTj ) we estimate 
eVpk ~ x*&bX. When the Pdi s mechanism dominates 
(kp£ < 10 3 ) and T To, a power-law relation is ob- 

1 /3 

tained: Vpk oc [k^iP/g 2 ] . Similar arguments lead 






6 



FIG. 4. (a) Electronic temperature T vs. bias Vb for different 
pump power values P. The inset shows the dominant cool¬ 
ing pathways for different chemical potential and pump power 
values [for To = OK, k?£ = 100]. (b) Current J vs. Vb for 
different kp£ values, and P m 10/iW//xm 2 . The bias poten¬ 
tial at which current peaks (Vpk) and the peak-to-valley ratio 
(PVR) are sensitive to the amounts of disorder (see inset). 
Results shown in this figure (except the inset of panel a) are 
obtained for the system undoped at Vb = 0, as in Fig [I] 

to a disorder-controlled peak-to-valley ratio (PVR). This 
behavior is illustrated in FigQb) inset. 

We note that the NDR features may be somewhat 
smeared out by statistical fluctuations induced by disor¬ 
der or inhomogeneities. However, we do not expect these 
effects to destroy NDR. Indeed, optical heating occurs in 
/i m-wide areas and a typical carrier density is 10 12 cm -2 . 
At the same time, charge inhomogeneity lengthscales in 
graphene/hBN are a few tens of nm, whereas typical den¬ 
sity fluctuations are as low as ~ 10 11 cm -2 [28l 129] . 

Summing up, vertically-stacked graphene heterostruc¬ 
tures afford a platform to realize and explore a range of 
interesting optoelectronic phenomena due to photogen¬ 
erated hot carriers. One such phenomenon is the light- 
induced NDR effect discussed above, manifesting itself 
through the I — V dependence, acquiring an N-shaped 
character under optical pumping. Vertical heterostruc¬ 
tures use the full graphene area as a photoactive region, 
and possess a large degree of tunability. These proper¬ 


ties make the NDR effect potentially useful for design¬ 
ing new types of optical switches and photodetectors. 
Our estimates show that the NDR regime, facilitated by 
graphene’s unique optical and thermal properties, can be 
readily accessed in wide-barrier heterostructures. 
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SUPPLEMENT: ELECTRONIC COOLING 
PATHWAYS 


Here we summarize the main results on electron-lattice 
cooling in graphene, following Refs. 0 0 [26]. Elec¬ 
tron cooling in graphene is usually assumed to be dom¬ 
inated by three main mechanisms: acoustic and opti¬ 
cal phonon emission m, and disorder-assisted acous¬ 
tic phonon emission (“supercollisions”) [26]. We use the 
two-temperature model describing the electron and lat¬ 
tice subsystems by two different temperatures, T and Tq. 
This model is valid when the electronic system is ther- 
malized quickly due to fast carrier scattering, whereas 
the electron-lattice cooling occurs on a longer time scale. 

The contribution to cooling power due to acoustic 
phonons, obtained for pristine graphene, is given by 0 

POO 

P.C = lac (T - To) / dvvH if(u) + 1 - /(-!/)] , (15) 
Jo 


where /(e) = [e^ 6- ^ + l] -1 and the quantity in the 
prefactor 7 ac = depends on the electron-phonon 

coupling strength 0 . Here D is the deformation poten¬ 
tial, p is the mass density of the graphene monolayer. 
An explicit dependence on chemical potential p can be 
obtained in the degenerate limit /3p 1. In this case, 

the integral in Eq. (15) yields 


7acF(T - T 0 ). 


(16) 


This contribution to cooling, due to its strong depen¬ 
dence on p, becomes very small near the Dirac point. 

Disorder-assisted acoustic-phonon cooling originates 
from electron-phonon scattering in the presence of dis¬ 
order, such that part of phonon momentum is absorbed 
by disorder. Evaluated for a short-range disorder model, 
this mechanism yields cooling power m 

o n 2 p 

ft.-w.v-u), < 17 > 

where kpl is the dimensionless disorder mean free path 
parameter, s is the speed of sound, and the degenerate 
limit f3p 1 is assumed. The quadratic dependence of 
Pdis on p means that this contribution can win over P ac 
near the Dirac point. 

Using values of R ^ 20eV, s ^ 2 • 10 4 m/s, p ~ 
7.6-10 -11 kg/cm 2 , we estimate 7 ac ~ 0.5 pW/pm 2 K(eV) 4 
and 7di S - 5 • 10“ 4 /iW/(eV) 2 K 3 /im 2 . 

Lastly, the cooling power for the optical phonon path¬ 
way equals 0 

Popt = 4^ H ^ 3 ) 2 [Nel{uJo) ~ N ph(wo)] HT ' m) ’ (18) 


where a flat optical phonon dispersion with Tiujq = 0.2 eV 
is assumed, and a = 1.42 A is the interatomic distance. 






The quantities N e \(uJo) and Aph(^o) represent the Bose 
distribution [e^ huJo — l] -1 evaluated at the electron and 
lattice temperature, T and To, respectively. The quantity 
T(T, /i) is a dimensionless integral 


/ oo 

dx\x(x - 1)| [f(hw 0 (x - 1)) - f(hu 0 x )], 

-OO 

(19) 

with /(e) defined above directly after Eq.(15). For weak 
doping, p < ftcjo, and k^T Tiujq, we can approximate 
/(e) by a step function 0(—e). Integration in Eq. (19) 
then yields T(T, /i) « 1/6, giving 


7-* TAT / N AT / M h(hUJo) 3 

P opt - 7opt [^Vel(^o) - ^Vph(^o)] , 7opt - 2A7Tpa A (hv) 2 ' 

( 20 ) 

At temperatures T, T 0 <C this gives a simple expo¬ 
nential dependence: 


Popt = 7opt [e-^ 0 - . (21) 

The exponential temperature dependence in Thpt makes 
this contribution small at low temperatures. However, 
the value P Qpt grows rapidly with tempertaure, becoming 


large at temperatures k^T > 0.27kj 0 . In addition, even 
at low temperatures, the effect of P Qpt can be important 
near the Dirac point, where other contributions are small 
since P ac oc //, P dis oc fi 2 . 

Using the above results we estimate the carrier densi¬ 
ties that mark the onset of the hot-carrier regime under 
typical experimental conditions. Reaching these carrier 
densities by adjusting by the gate potential and voltage 
bias triggers the NDR regime (see main text). Below we 
assume, for concreteness, that electron-lattice cooling is 
dominated by the disorder-assisted mechanism. Taking a 
typical pump power value P ~ 10/iW//im 2 , equating it 
to T d is, and using values kpl ~ 100 and T = 2T 0 ~ 600 K, 
Eq. ( [T7| ) yields typical carrier densities on the order of 
n ~ l(r 2 cm -2 . Similar values are obtained for the acous¬ 
tic phonon mechanism in the absence of disorder. The 
hot-carrier transport regime is therefore realized when 
the carrier density is n < 10 12 cm -2 , whereas field emis¬ 
sion dominates for n > 10 12 cm -2 . Such doping values 
are routinely obtained in graphene systems through elec¬ 
trostatic doping, indicating that the densities required 
for NDR are easily reachable under realistic experimen¬ 
tal conditions. 





